Temperature Anisotropy in a Driven Granular Gas 
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When smooth granular material is fluidized by vertically shaking a container, we find that the 
temperature in the direction of energy input always exceeds the temperature in the other directions. 
An analytical model is presented which shows how the anisotropy can be traced back to the inelas- 
ticity of the interparticle collisions and the collisions with the wall. The model compares very well 
with molecular dynamics simulations. It is concluded that any non-isotropic driving of a granular 
gas in a steady state necessarily causes anisotropy of the granular temperature. 
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As compared to an ordinary, molecular gas, the hall- 
mark of a granular gas is its permanent dissipation of 
energy due to inelastic collisions. Whereas an isolated 
molecular gas will sustain its motion for an infinite 
amount of time, the only true equilibrium state of gran- 
ular matter is the one where it is at rest. Hence a steady 
supply of energy is required to keep a granular gas alive, 
giving rise to prototypical non-equilibrium systems with 
many striking phenomena [1, 2]. The one addressed in 
this Letter is the crucial temperature anisotropy within 
a granular gas. It is observed to be a significant effect 
in both numerical simulations [3-8] and experiments [9- 
12]. Although it has been studied in the context of a 
random restitution coefficient model [13, 14], a theoreti- 
cal explanation is still lacking. Here we provide such an 
explanation, in which for analytical convenience we will 
restrict ourselves to a dilute granular gas, fluidized by 
vertically vibrating a container. 

So what causes the anisotropy? We approach this ques- 
tion by a theoretical model in combination with event 
driven molecular dynamics (MD) simulations, and show 
that the anisotropy results from the following characteris- 
tic feature of such a gas: The distribution of energy from 
the vibrating bottom towards the horizontal directions 
occurs through the very same mechanism that also con- 
stitutes one of the major sources of energy dissipation, 
i.e., the collisions between the particles. This result car- 
ries over to any granular gas with a non-isotropic energy 
source. 

The setup we will consider in our present work consists 
of a granular gas in a container with a square-shaped bot- 
tom of side length L in the x-y-planc and infinitely high, 
vertical side-walls. Gravity acts with g = 9.81 m/s 2 and 
the gas is fluidized by vertical vibrations of the bottom 
about z — with amplitude a and frequency / - typically 
of triangular (piecewise linear, symmetric) or sinusoidal 
shape. The gas consists of N identical hard spheres with 
radius r p and mass m. We restrict ourselves to the case 
[3, 15-20] that only normal restitution [1, 2] contributes 
to the dissipative processes, with restitution coefficients 
e p for particle-particle and e w for particle- wall collisions, 
while collisions with the vibrating bottom are taken to 
be perfectly elastic. 

After initial transients have died out, we expect (and 
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FIG. 1: Solid lines: Temperature components (1) from MD 
simulations of the driven granular gas model as specified in the 
main text with sinusoidal driving and parameters / = 72.648 
Hz, a = 1.00 mm, TV = 50, r p = 1.18 mm, and L = 2.265 
cm. (a): e p — 0.9 and e w — 1 (reflective sidewalls). (b): 
e p = 0.95 and e w = 0.95 (dissipative sidewalls). Dashed lines: 
Counting only particles that move downward, i.e. restricting 
the average in (1) to negative u z -values. At large z the statis- 
tics deteriorates since particles rarely reach such heights. The 
curves T x (z) and T y (z) practically coincide. 



observe) a stationary probability distribution p(r, v) of 
the particle positions r := (x,y,z) and velocities v := 
(v x ,v y ,v z ). The quantities of central interest are the 
temperature components 
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where i £ {x,y, z} and k is Boltzmann's constant. This 
Ti(z) is directly proportional to the average kinetic en- 
ergy of the particles in a horizontal layer at height z in 
either of the three spatial directions i. 
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In evaluating our MD simulations we replace ensem- 
ble averages in (1) by time averages (justified by ergod- 
icity) and work in units with m = 1 and k = 1. A 
representative result is depicted in Fig. 1: As expected 
for symmetry reasons, the horizontal temperature com- 
ponents T x (z) and T y (z) are practically indistinguish- 
able. In contrast, the vertical temperature component 
T z (z) is significantly larger than T x (z) and T y (z). For 
perfectly elastic particle-wall collisions (Fig. la) the z- 
dependcnce of the temperature components is very weak 
except for a region directly above the bottom. There, 
the energy-input by the driving leads to increased up- 
ward particle velocities, as shown by the dashed lines in 
(Fig. la). For inelastic particle- wall collisions (Fig. lb) 
the z-dependence of the temperature components is more 
pronounced. Yet, in both cases the differences between 
vertical and horizontal temperatures are more important 
than the z-dcpcndcnccs. 

Our theoretical analysis of the observed temperature 
differences starts with the well-known conservation laws 
of energy and momentum for a dilute granular gas, de- 
rived from Boltzmann's equation [15-17]. For a station- 
ary system, they read, in terms of the local heat flux J(r) 
and stress (or pressure) tensor P(r): 



V • J(r) = /(r) , V • P(r) = n(r)f (r) 
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Here, I(r) is the local energy dissipation rate per unit 
volume, n(r) is the local particle density, and the force 
f(r) = —mge z . Integrating the first equation (2) over 
the container volume V we obtain with Gauss' theorem 
that the energy dissipation rate due to particle-particle 
collisions Q pp := J v I(r)dr must be equal to the total 
flux of energy through the boundaries. The latter can 
be decomposed in the influx Q in := L 2 J z (0) of energy 
through the vibrating bottom of area L 2 and the energy 
dissipation rate Q w due to particle-wall collisions. This 
gives 
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Crucial to the present model is that the temperature 
components defined in (1) are treated separately. To 
our knowledge, all existing theories for driven granular 
gases in a steady state without net flow of material are 
based upon the assumption of an isotropic temperature, 
and many of them also on isotropic stress [15-18, 28]. 

The temperature components T, are related to the di- 
agonal elements of the stress tensor by a generalized ideal 
gas law Pa = nkTi [15-17]. Motivated by our MD simu- 
lations, we assume that each temperature component Ti 
is approximately constant within the entire container vol- 
ume (see also [19, 21]). Because of symmetry the stress 
tensor is diagonal and T x = T y =: 2^., with which the 
second equation (2) can now be readily integrated to yield 
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Furthermore, as exemplified by Fig. 2, our MD simula- 
tions show that the particle velocity components can be 
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FIG. 2: Horizontal (solid red and blue dots) and vertical 
(open, black dots) velocity distributions of particles in a layer 
between 2 = and z — 5 mm obtained by MD simulations 
for the same system as in Fig. la but with / = 30 Hz and 
L = 4.8 cm. The solid lines are Gaussian fits [29]. 



assumed as Gaussian distributed in very good approxi- 
mation for a wide range of parameter values [22, 29]. All 
together, we thus arrive at the following approximative 
distribution function for particle position and velocity: 
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The next main idea is to determine the two unknowns 
T z and T^, in (5) by means of two energy balance rela- 
tions. The first of them is (3). To obtain the second, 
we observe that particle-particle collisions not only cause 
energy dissipation but also a transfer of kinetic energy 
from the horizontal direction into the vertical direction, 
and vice versa. In the steady state the net effect must 
be an average loss of kinetic energy per time unit Q z in 
the vertical direction which is exactly balanced by the 
incoming energy flux Qi n through the vibrating bottom: 



Qin — Q z 



(6) 



The remaining, rather technical task is to explicitly 
determine all the energy fluxes Q appearing in (3) and 
(6) with the help of the approximation (5). In order to 
evaluate Q z , we first note that the change of kinetic en- 
ergy in the vertical direction in a single particle-particle 
collision is 
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where Vj and v1j are the velocities of the two colliding 
particles (j = 1,2) before and after the collision, respec- 
tively. Due to our assumption that only normal resti- 
tution contributes to the dissipative processes, we have 
v* = Vj + (— 1) J (1 + e p )[(vi — v 2 ) ■ n]n/2, where n is 
the collision normal vector. To determine Q z one es- 
sentially has to introduce this result for v* into (7) and 
then average according to (5). More precisely, first q z in 
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(7) is multiplied by the collisional volume per unit time 
(i/2)7r(2r p ) 2 |i>2 — «i|5(|t"i — r 2 | — 2r p ), where the factor 
1/2 arises since collisions only can happen if (vi— v\)-n < 
0. Next we multiply with (l/2)p(ri, Vi)p(r2, V2) accord- 
ing to (5) and integrate over all Vj and rj (within the 
container volume). The factor 1/2 is needed since every 
collision appears twice in the above considerations. This 
gives, after a substantial amount of algebra: 



(1 + e p ) (Nr p /L) 2 g^/irmkT z F(T har /T z ) (8) 
(8-4e p )s- (l + e p )t 



ds / dt y/s + t 
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A similar averaging of the total energy loss in a single 



particle-particle collision q pp 
yields 
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In the same spirit one can evaluate the total dissipation 
rate due to particle-wall collisions Q w with the result 

Q w = A{\~e 2 w )N{kT lva ,f/ 2 /{L^m) (10) 

Finally, a somewhat lengthy but straightforward calcula- 
tion yields the following expression for the energy input 
rate at the perfectly elastic, vibrating bottom of the con- 
tainer: 
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where 7 := y / 2kT z /mu 2 , u =: n\/2af for sinusoidal 
driving of the bottom, and u := Aaf for triangular driv- 
ing. In both cases, for large 7 (which is the typical situa- 
tion in the dilute systems under study) the function ^(7) 
approaches unity [30]. 

In the absence of wall dissipation (e w = 1), Eqs. (3), 
(6), (8)-(ll) imply 



FiT^/T,) = (1 - e p )G(T^/T z ) , [e v 



(12) 



Closer inspection shows that for any < e p < 1 a unique 
solution < TiKf/Tz < 1 of (12) exists. In particular, for 



small 1 — e p one finds the leading order asymptotics 



WT 2 = 1 - (5/2) (1 - e p ) , [e t 



(13) 



Thus for e p < 1 and perfectly reflecting walls, the model 
predicts that the horizontal temperature is always 
smaller than the vertical temperature T z . Moreover, the 
ratio Tha- /T z solely depends on e p but not on any details 
of driving shape and strength, particle density, particle 
size, or compartment geometry. The comparison with 
MD simulations in Fig. 3 is excellent. The inset of Fig. 3 
shows that the remaining deviations can largely be at- 
tributed to finite size effects. 



FIG. 3: Ratio of horizontal and vertical temperatures ver- 
sus particle-particle restitution coefficient. Solid line: Theo- 
retical prediction (12). Dashed line: Theoretical asymptotic 
behavior (13). Symbols: MD simulations for the same sys- 
tem as in Fig. la with N — 20 (*) and N = 50 (•) parti- 
cles. For e p 7^ 0.9, a different choice for the frequency / such 
that the T z was (approximately) the same for all e v . Inset: 
Same, but for fixed e p = 0.9 and variable local solid fraction 
a := (4/3)7n-p max[n(r)], either by varying particle number 
N (0, at r = 1.18 mm) or radius r p (□, at N = 50). 



In the general case one obtains two transcendental, al- 
gebraic equations for the two unknowns and T z by in- 
troducing (8)-(ll) into (3) and (6). While existence and 
uniqueness of solutions can still be demonstrated analyt- 
ically, their quantitative determination is only possible 
numerically. An example is depicted in Fig. 4, compar- 
ing very well with MD simulations. As expected, dissi- 
pative walls tend to reduce the horizontal temperature 
Tfror since they add a source of dissipation for the hori- 
zontal kinetic energy. If we increase the particle num- 
ber, the gas becomes more dense, and the number of 
particle-particle collisions will increase much faster than 
the number of particle- wall collisions. Therefore, will 
first increase, in sharp contrast to the overall tempera- 
ture [T z + 22/JQTJ/3 which must decrease with increasing 
particle density. Eventually, particle-particle collisions 
will dominate the system, and /T z will asymptotically 
tend to the value it would have with reflecting walls. 

In conclusion, we have numerically observed large dif- 
ferences between the vertical and horizontal tempera- 
tures in vertically driven granular gases subjected to 
gravity. We introduced a theoretical model based on 
an approximative Maxwcll-Boltzmann distribution with 
anisotropic but homogeneous temperature (5), justified 
by our MD simulations. Both for reflecting and dissipa- 
tive walls of the container we find that the theoretical 
model gives good quantitative agreement with the simu- 
lations. 

The difference of the horizontal and vertical temper- 
atures from the isotropic value is a significant correc- 
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FIG. 4: Horizontal (blue) and vertical (red) temperatures ver- 
sus particle number. Lines: Theoretical predictions from (3), 
(6), (8)-(ll) for triangular driving and parameters / = 30 
Hz, a = 1.00 mm, r p = 1.18 mm, L = 4.8 cm, e p — 0.95 and 
e w — 0.95. Symbols: MD simulations. Inset: Same data in a 
log- log plot. 

tion, at least of the same order as those resulting from 
taking into account the non-constancy of the tempera- 
ture and density profiles, or from embedding Chapman- 
Enskog corrections to the velocity distributions into the 
theoretical framework [15, 16]. 



At the root of the temperature anisotropy lies the fact 
that the transfer of kinetic energy between different spa- 
tial directions and its dissipation arise out of the same 
mechanism: the collisions between the particles. Thus, 
the anisotropy of the temperature is a necessary conse- 
quence of the anisotropy of the driving. The present 
work indicates that one may obtain an improved hydro- 
dynamic description for dilute granular gases by starting 
out from an anisotropic velocity distribution (instead of 
an isotropic one) in deriving hydrodynamics equations 
from Boltzmann's equation. 

The basic concept of the model applies to any situation 
in which the energy-input in a granular gas is anisotropic, 
always predicting a higher kinetic energy content in the 
main direction of energy- input. 
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